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' The adiabatic self-consistent collective coordinate (ASCC) method is a practical micro- 

scopic theory of large-amplitude collective motions in nuclei with superfluidity. We show 
that its basic equations are invariant against transformations involving the gauge angle in 
the particle-number space. By virtue of this invariance, a clean separation between the large- 
■^j- ' amplitude collective motion and the pairing rotational motion can be achieved, enabling us 

to restore the particle-number symmetry broken by the Hartree-Fock-Bogoliubov (HFB) 
approximation. We formulate the ASCC method explicitly in a gauge-invariant form. In 
solving the ASCC equations, it is necessary to fix the gauge. Applying this new formulation 
to the multi-0(4) model, we compare different gauge-fixing procedures and demonstrate that 
£ — , calculations using different gauges indeed yield the same results for gauge-invariant quanti- 

ties, such as the collective path and quantum spectra. We suggest a gauge-fixing prescription 
that seems most convenient in realistic calculations. 
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§1. Introduction 

Construction of the microscopic theory of large amplitude collective motion is a 
long-standing and fundamental subject of nuclear many-body problem. 1 )~ 35 ) As is 
well known, pairing correlations play crucial roles in low-lying states of medium and 
heavy mass nuclei, and they are taken into account in the Hartree-Fock-Bogoliubov 
(HFB) mean-field theory by breaking the particle number conservation. 36 )~ 38 ) The 
broken particle-number symmetry can be restored, however, by making the self- 
consistent quasiparticle random-phase approximation (QRPA) in addition to the 
HFB mean field. 36 )" 38 ) It is an advantage of the QRPA that number-conserving col- 
lective modes, such as shape-vibrational modes, are exactly decoupled from number- 
fluctuation modes. The latter modes are associated with the nucleon- number degrees 
of freedom and called the pairing rotational modes. It is a unique feature of nuclei as 
finite quantum systems that such a rotational motion in the gauge space is actually 
observed in quantum spectra. 39 ) Since the applicability of the QRPA is limited to 
small-amplitude collective motions, it is highly desirable to extend the QRPA to a 
general theory keeping its decoupling feature. Such a theory should be capable of 
describing the interplay between large-amplitude collective motions and the pairing 
rotational modes. 

The self-consistent collective coordinate (SCC) method 7 ) is a microscopic the- 
ory of large-amplitude collective motion based on the time-dependent Hartree-Fock 
(TDHF) method, which enables us to extract the collective submanifold in a fully 
self-consistent manner. The SCC method was originally formulated 7 ) for systems 
without pairing correlations, and then extended 18 ) to systems with superfluidity. To 
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extract the collective submanifold embedded in the time-dependent Hartree-Fock- 
Bogoliubov (TDHFB) phase space, the number and angle variables describing the 
pairing rotational motion are explicitly introduced 18 - 1 in addition to the collective 
variables describing the large-amplitude collective motion. This extended version 
of the SCC method has been applied successfully to various kinds of anharmonic 
vibration and high-spin rotational motion. 40 ^ 51 -' However, its solution relied on an 
expansion technique with respect to the collective coordinates and momenta around 
the HFB states. Thus, it is difficult to describe collective motions with a genuine 
large-amplitude nature. Recently, the adiabatic SCC (ASCC) method has been pro- 
posed to overcome this difficulty 52 ) The ASCC method is an alternative way to 
solve the SCC basic equations assuming that the large-amplitude collective motion 
of interest is slow (adiabatic). Under this assumption, the basic equations of the 
SCC method are expanded up to the second order in the collective momentum, but 
no expansion is made with respect to the collective coordinate. This has been ap- 
plied 53 ) to shape-coexistence phenomena, 54 )~ 56 ) where the large-amplitude collective 
motions take place between the oblate and prolate HFB equilibrium shapes. How- 
ever, the calculation often leads to a numerical instability, caused by a redundant 
degree of freedom in the ASCC equations, which was treated with an additional 
condition by hand. 57 ) Now we understand the origin of this redundancy. It turns 
out to be due to a gauge invariance of the ASCC equations. 

The main purpose of this paper is to formulate the ASCC method in a way 
to manifest the invariance against a transformation with respect to the angle vari- 
able in the gauge space. This invariance is crucial to achieve a clean separation 
between the large-amplitude collective motions and the pairing rotational motions, 
and to restore the particle-number symmetry broken by the HFB approximation. 
This new formulation of the ASCC method also provides a justification to the pre- 
scription adopted in our previous work. 57 ) We apply the method to the multi-0(4) 
model 57 ) -62 ) with different choices of the gauge and test the internal consistency of 
the proposed scheme by carrying out a detailed numerical calculation. 

This paper is organized as follows: In £j2]the basic equations of the ASCC method 
are recapitulated. In $3] the gauge-invariant formulation of the ASCC method is 
given. It is applied to the multi-0(4) model in ^5] and tested by numerical calcula- 
tions in Sj5j Concluding remarks are given in $6) 

§2. Basic equations of the ASCC method 

2.1. Basic ideas 

Time evolution of large amplitude collective motion is described by the time- 
dependent variational principle, 

5(m\i§- t -Hm)) =0, (2-1) 

for the time-dependent HFB state vector \4>(t)). Assuming that the time-dependence 
of the collective motion is governed by the collective coordinate q and the momentum 
p, we parameterize the time-dependent HFB state vector as 

|0(*)> = \<P(q,p,<p,n)). (2-2) 

Here, (p represents the gauge angle conjugate to the particle-number n = N — Nq. 
We measure the particle number from a reference value Nq specified below, and 
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assume, for simplicity, that the number of collective coordinate is one. We define 
the intrinsic state vector \(j>(q,p,n)) in the particle-number (gauge) space by 

\<l>(q,p,<p,n))=e-*'*\if>(q,p,n)), (2-3) 

where N = N — Nq. Two sets of collective variables (q,p) and (<p,n) are determined 
such that the canonical variable conditions, 

d 

(<p(q,p,n)\i— \(j)(q,p,n)) 

(cj)(q,p,n)\ N \(j)(q,p,n)) 

are satisfied. Here the generating function S is an arbitrary function of q,p, (p, and n. 
We choose S = because it is appropriate to the adiabatic approximation. 15 )' 16 )' 20 ) 
The collective Hamiltonian is defined by 

H(q,p,n) = (4>(q,p,ip,n)\H\4>(q,p,ip,n)} = (4>(q,p, n) \ H \ <ft(q,p, n)) . (2-5) 

Note that it is independent of the gauge angle tp, because the Hamiltonian commutes 
with the particle-number operator N. 

The equation of collective path is obtained by replacing the time derivative term 
in Eq. (|2Tp with derivatives with respect to four collective variables 

6 «,, P ,„)| H - , - __ + --Nj W ,,p,„)) = 0, (2-6) 

where the canonical equations of motion for collective variables (q,p) and (ip, n) are 
used in order to eliminate the time-derivative of the collective variables. Equations 
(E31), (I23D, and ([M]) constitute the basic equations of the SCC method. 7 )' 18 ) 

2.2. Basic equations of the A SCC method 

Assuming that the large-amplitude collective motion is slow, i.e., p is small, let 
us write the TDHFB state vector \<f>(q,p, n)) in the following form: 

\0(q,p,n))=e^ +m °^ \<f>(q)) . (2-7) 
Here Q(q) and 0(q) are Hermitian one-body operators, which may be written as 

Q(<l) = X] {Q^(l) a i a l + Qafi{q)*apac^ , (2-8) 

a/3 

0{q) = i (e«0te)44 " ^(Q)*apa a ) , (2-9) 

af3 

where aj, and a a are quasiparticle creation and annihilation operators associated 
with the time-even state vector \(p(q)} and satisfy a a \4>(q)) = 0, and n = N — N , 
Nq being the expectation value of the particle number with respect to \4>(q)). We 
shall discuss in ^}3]that it is also possible to adopt a slightly different representation 
for Q(q). 

Substituting (|2-7p into (|2 • 4[) and comparing the coefficients of the zero-th and 
the first order terms in p and n, we obtain the canonical variable conditions in the 



dS a OS 

P + -x~> (H ( liP, n )\ t^- \<t>(q,P,n)) = -— , 
oq top op 



dS d 



dS_ 

dn ' 
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adiabatic limit: 



(0(g) | P(g) 

M?)IQ(?) 

(0(g) | iV 
(0(g)|0(g) 
(0(g)|[@(g),iV] 
(0(g)|[Q(g),0(g)] 



0(g)) = o, 
0(g)) = o, 
0(g)) = o, 
0(g)) = o, 
0(g)) = h 
0(g)) = o, 

0(g)) = -] 



where -P(g) is the local shift operator defined by 



(2-10) 
(2-11) 
(2-12) 
(2-13) 
(2-14) 
(2-15) 

(2-16) 



(2-17) 



Differentiating (|2- 1 1|) and f)2- 12j) with respect to q and using (|2- 16|) . we obtain 

(0(g)|[Q(g),%)] 10(g)) =*, (2-18) 
(0(g)| [N,P(q)} |0(g)) =0. (2-19) 

The collective Hamiltonian f)2 • 5[) is also expanded up to the second order in p 
and the first order in n, 



H(q, p, n) =V{q) + -B{q)p 2 + \(q)n, 



(2-20) 



where the collective potential V(g), the inverse mass parameter B(q), and the chem- 
ical potential A(g) are defined by 



V(q) =H(q,p,n) 
ld 2 H 



p=n=0 



(0(g)|il |0(g)) 



B{q) 



2 dp 2 



p=n=0 



(4>{q)\[[H,iQ{q)],iQ{q)] |0(g)) 



p=n=0 



(0(g)|[il,*<9(g)]|0(g)). 



(2-21) 
(2-22) 
(2-23) 



We obtain the ASCC equations by expanding the equation of collective path 
(!2-6h with respect to p and n, and requiring that the variations vanish for each order 
in p and n. In the zero-th order, we obtain the moving-frame HFB equation: 



where 



<5(0(g)|tf M (g)|0(g))=O, 



H M (q)=H-\(q)N-^-Q(q), 



(2-24) 



(2-25) 



is the moving-frame Hamiltonian. In the first and the second orders, we obtain the 
local harmonic equations (also called the moving- frame QRPA equation): 



5 (0( g )| [H M (q),iQ(q)] - B{q)P{q) |0(g)) = 0, 



(2-26) 
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8(t(q)\[H M (q),P(q)]-iC(q)Q(q) 



1 ptW f)\ 

[[H M (q), -5-Q(g)] s »Q(g)] " *^ l<K<?)> = 0, (2-27) 



where 

„, , <9 2 F 1 dBdV , 

Note that in Ref. 52) the curvature term l/2B(q)[[H M (q),dV/dqQ(q)],iQ(q)\ is 
linear lized with respect to Q(q) using the relation 

BV 

(H - X(q)N) A = —Q(q), (2-29) 

where the superscript A in Eq. (|2-29|) denotes the two-quasiparticle creation (a' a') 
and annihilation (aa) part of the operator in the parenthesis. Hereafter, we call 
them "^4-part" and the a^a terms "B-part". The collective variables (q,p) and the 
collective Hamiltonian 7i(q,p,n) are determined by solving the ASCC equations, 
(|2-24p . (|2-26p . and (|2-2T[) . under the canonical variable conditions. Note that we can 
make a scale transformation of the collective coordinate q such that B(q) = 1. We 
adopt this choice. Then, C(q) represents the curvature of the collective potential: 

d 2 V(q) 

C(q) = (2-30) 

§3. Gauge invariance of the ASCC equations with respect to the 
pairing rotational degree of freedom 

3.1. Gauge invariance at the HFB equilibrium point 

As mentioned in the preceding section, the first step to solve the ASCC equations 
is to find a solution at one of the HFB equilibrium point denoted by q = qo, which 
corresponds to the local minimum of the collective potential V{q) satisfying dV/dq = 
0. The moving-frame HFB equation reduces to the conventional HFB equation at 
the equilibrium point, 

8(4(q )\H-\(q )N\<P(q ))=0. (3-1) 
The local harmonic equations at the equilibrium point are given by 

8 {<f>(q )\ [H - X(q )N, iQ(q )] - B(q )P(q ) |0(g o )) = 0, (3-2) 
5 (<f>(go) | [H - X(q )N, P(q )} - iC(q )Q(q ) - i—N \<j>(q )) = 0. (3-3) 

These equations reduce to the QRPA equations if the quantity dX/dq vanishes. In 
other words, the QRPA solution corresponds to the special solution of the local 
harmonic equations with dX/dq = 0. 

Let us consider the following transformations: 

Q(qo)^ Q(qo) + aN A (q ), (3-4a) 
0(qo) - ©(go) + aP(q ) (34b) 

^(go) ^(go) ~ aC (qo)- (3-4c) 
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Here, a is an arbitrary number and N A denotes the ^4-part of the number operator 
N. This is a kind of gauge transformation with respect to the pairing rotational 
degree of freedom. We can easily confirm that the local harmonic equations at the 
HFB equilibrium point, (|3-2p - (|3-3p . and the canonical variable conditions, f)2- 10j) - 
(|2-15p and (|2- 18[) - (|2- 19[) . are invariant under this transformation. Due to this invari- 
ance, the solution of the local harmonic equations is not uniquely determined at the 
HFB equilibrium point. If we choose a value of a such that d\/dq = holds, the 
local harmonic equations coincide with the conventional QRPA equations. We can 
choose other values of a, however, if it is more convenient. 

3.2. Gauge invariance at non- equilibrium points 

At non-equilibrium points, dV/dq is non-zero and the moving- frame Hamilto- 
nian (|2-25|) depends on the collective coordinate operator Q(q). Still we can general- 
ize the above consideration at the HFB equilibrium point to a general off-equilibrium 
point q on the collective path: It is straightforward to confirm that all the basic 
equations of the ASCC method [i.e., the collective Hamiltonian TC(q,p,n), (|2-20j) . 
the inverse mass parameter B(q), (|2-22|) . the moving-frame HFB equation, (|2-24j) . 
the local harmonic equations, (|2-26p - (|2-27p . and the canonical variable conditions, 
(|2-10p - (|2-15p and (|2-18p - (j2-19p ] are invariant under the following transformations 
with respect to the pairing rotational degree of freedom: 

Q{q) Q{q) + aN, (3-5a) 
0(q)^0(q) + aP(q), (3-5b) 

A(g)- X(q)-a^(q), (3-5c) 
^(q)^^(q)-aC(q), (3-5d) 

if the collective coordinate operator Q(q) is constructed such that it is exactly com- 
mutable with the number operator N, 

[Q(q),N} = 0. (3-6) 

In association with the above transformations of Q(q) and 0(q), the original 
TDHFB state vector 

\<Kq,p,<p,n)) = e -^e^e^ \cf>(q)) , (3-7) 

is transformed as 

\(/){q,p,ip,n)) ^ e -ivN e ip(Q( q )+aN) ei n(0( q ) +a P( q )) 

=e -i( V -ap)N e ipQ(q) e ine(q) ^ g _ m ^ _ (3.3) 

Here, the relation 

\Hq + Sq)) = e- i5qp(q) \<t>(q)) (3-9) 

is used and the operators, 0{q) and P(q), are treated as commutable with each other 
under the adiabatic approximation. We also note that the expression (|3 7[) is slightly 
different from (|2T|) : the difference between e ipQ(g)+in0(g) and jvQ{q) e ine{q) ^ however, 
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gives rise to only higher order contributions, which are ignored in the adiabatic 
approximation under consideration. We see that the gauge angle ip changes to ip— ap 
due to the transformation f)3- 5[) . Thus, hereafter, we briefly call the transformations 
(|3 5j) "gauge transformations" , independence of the choice of a "gauge invariance" , 
and choice of the value of a "gauge fixing." 

The commutability, (|3-6p . implies that Q(q) is a normal one-body operator writ- 
ten in terms of the nucleon creation and annihilation operators (c^ , c) in the following 
form: 

Q(q) =Y t QiM ■ c\c, : = Q A (q) + Q B (q). (3-10) 

Here, the symbol : : denotes the normal product part when Q(q) is written in terms 
of the quasiparticle operators (a', a) defined at \<fi(q)), and the coefficients satisfy 
the relation Qij(q) = Qji(q)*, because Q(q) is supposed to be Hermitian. Note that 
this Q(q) operator contains the 5-part, Q B (q) = Ylki Qki a \ a ii as wen as ^ ne 
Q A (o) = J2m Qkl a k a i nx ' Accordingly, the relation (|2-29p does not hold for this 
Q(q) operator. 

In this way, we arrive at a new formulation of the ASCC method in which the 
gauge-invariance (|3-5p is manifest. The gauge-invariant ASCC method consists of 
the basic equations which are the same as those in the orig inal ASCC method 52 ) 
except for the use of Eq. (|2-27p and the Q(q) operator given in the form of (|3T0p . 

3.3. Gauge fixing and numerical algorithm 

The fact that the ASCC equations are invariant against the gauge transfor- 
mations ()3-5p indicates the necessity of choosing a particular gauge for numerical 
computation: If the gauge is not fixed, an instability with respect to the gauge 
degree of freedom might occur during the course of numerical calculation. Let us 
outline the procedure of gauge fixing and numerical algorithm for solving the gauge 
invariant ASCC equations. We start the calculation by solving the moving-frame 
QRPA at one of the HFB equilibrium points. A solution of the moving-frame QRPA 
at the HFB equilibrium point can be obtained, as discussed in $\3.1\ by choosing the 
gauge dX/dq = 0. Hereafter we call this gauge "QRPA gauge", because under 
this gauge, the moving-frame QRPA equations at the HFB equilibrium reduce to 
the conventional QRPA equations. As we shall see later, the numerical calculation 
using the QRPA gauge encounters a difficulty at inflection points of the collective 
potential V{q). It is possible, however, to choose another gauge that is free from 
this difficulty. With the use of the multi-0(4) model, we shall explicitly show in §5\ 
how this is done. 

Since the moving- frame HFB equation at non-equilibrium points contains Q(q) 
that should be determined by the local harmonic equations, we have to resort to 
an iterative procedure. We proceed to the direction of the lowest energy solution 
of the moving-frame QRPA, and successively find the solutions in the following 
manner. Suppose that we already obtain the solution at q' = q — 5q where 5q is the 
numerical mesh size in computation. The moving-frame HFB equation at q for the 
n-th iteration 

^(?))=0, (3-11) 
H$\q) = H — X^(q)N - ^ (q)Q^)(q), (3-12) 



5U n \q) B®(q) 
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is solved under the following two constraints 



(q) 



N 



(q) 



An 



4> {n \q) Qiq-Sq) ^\q)\ = Sq, 



(3-13) 
(3-14) 



which are derived from the canonical variable conditions, (|2- 12j) and (|2-16p . respec- 
tively. In starting this iterative procedure at q, the neighboring solution Q(q — Sq) 
(or a linear combination of the moving-frame QRPA modes at q — 5q) may be used 
as an initial trial for the operator Q(q)^. The moving-frame QRPA equations for 
the n-th iteration are written as 



5U n \q) [H$(q),iQ (n Hq)} ~ B^(q)P (n) (q) = 0, 



(3-15) 



S(^ n \q)\ [Hi}\q),P {n) (q)] ~ ^[tf(g),^(g)],#)] 

f) \ (n) . 
- iC^(q)Q^(q) - i— (q)N ^ n \q)) = 0. (3-16) 

As the curvature term is nonlinear with respect to Q(q), we replace one of the 
operator Q(q) with that of the previous iteration step Q^^iq)- This procedure 
will be discussed in detail for the multi-0(4) model in £14.51 Thus the moving- frame 
QRPA equations are linearlized with respect to Q^ n \q) and p( n \q). Its A-part, 
Q A ( n \q), can be determined in a manner similar to solving the conventional QRPA 
equations. After obtaining Q A ^ n \q), its -B-part, Q B ^ n \q), is determined with the 
requirement, Eq. (|3-6p . This iterative procedure is repeated until we get convergence 
for X {n) (q),Q < - n) (q), and |0 (n) (g)>- 

§4. Application of the gauge-invariant ASCC method to the multi-0(4) 

model 

In the following sections, we apply the scheme formulated above to the multi- 
0(4) model, and discuss gauge-fixing conditions suitable for solving the gauge- 
invariant ASCC equations. We discuss excitation spectra and transition strengths in 
systems with definite particle number N = Nq, so that we put n = in the following 
sections. 

4.1. The multi-0(4) model 

The multi-0(4) model Hamiltonian has been used to test the validity of mi- 
croscopic theories of nuclear collective motion. 57 ^ 62 ) We employ exactly the same 
model Hamiltonian as was used in the previous work. 57 ) The model is constructed 
in terms of the generators of the 0(4) symmetry, 



A ) = E 4™ C i-™> B j = E a irnc) m c)_ m , (4-1) 

m>0 m>0 

Nj = c im c 3mi Dj = ^ a jm c -j m c jmt (4- 2) 
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where the nucleon creation and annihilation operators (cj m ,Cj m ) are used and the 
quantity Oj m is given by 

ajm \ -1 \m\ > [ij/2. 1 > 

These four operators represent the monopole pair, the (simplified) quadrupole pair, 
the particle number, and the (simplified) quadrupole operators for each j-shell, 
respectively. The model Hamiltonian is written in the following form: 

H = h - ^G (AU + iit) - ^G 2 (&B + B&) - \ X D\ (4-4) 

j 

where defined by 

Ai = J2 A P ^ = J2djB], N = J2N h D^djDj, (4-5) 
j j j j 

and dj represents the quadrupole matrix element. The first term on the right- 
hand side of Eq. (|4-4p is the single-particle Hamiltonian, giving a spherical single- 
particle energy for each j-shell which possesses (2i7j)-fold degeneracy (2i?j = 
2j + l). The other terms represent the residual two-body interactions: the monopole- 
pairing interaction, the quadrupole- type pairing interaction, and the quadrupole- 
type particle- hole interaction. Their interaction strengths are denoted by Go, G*2, 
and Xi respectively. Here, the operators A and N are the monopole-pair and the 
number operators, while B and D represent the simplified quadrupole-pair and 
quadrupole particle-hole operators, respectively. 

The residual interactions of this model are separable and we can always write 
such a separable Hamiltonian in the following form: 

H = ho - ~ £ + \ E *sF^Ft\ (4-6) 

s s 

where 

Fi ± )^(F s ±F])/2 = ±Fi ± )t. (4-7) 

The superscripts (±) indicate the Hermitian or anti-Hermitian nature of the bilinear 
operator F. The multi-0(4) model Hamiltonian under consideration contains three 
kinds of residual interactions. The indices s =1, 2, and 3 on the operators F and 
the interaction strengths k s indicate the monopole-pairing, the quadrupole-pairing 
and the quadrupole particle-hole interactions, respectively: F s= i = A, F s= 2 = B, 
F s=3 = D, Ki = 2G , k 2 = 2G 2 , and k 3 = x- 

4.2. Quasiparticle representation 

To solve the ASCC equations, it is convenient to use the quasiparticle basis 
locally defined with respect to the state \4>(q)) on the collective path. For the multi- 
0(4) model, the Bogoliubov transformation to the quasiparticle creation and anni- 
hilation operators, a\{q) and Oj(g), satisfying the vacuum condition, aj(g) \4>{q)) = 0, 
is written 



4(q) \ _ (ui{q) -Vi(q)^ ' -- T 



a-i(q) J \Vi(q) Ui(q) J \ c_ 



c 



(4-8) 
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Here, the indices ±i represent the set of angular momentum quantum numbers 
(j,±m). The Bogoliubov transformation is locally determined on the collective 
path by the moving-frame HFB equation for a given collective coordinate operator 
Q{q). " 

Using the quasiparticle bilinear operators 



Al( q ) = a]( q )aUq), 

N i(q) = 4(q) a i(q) + aii(?)a-*(9)> 

the nucleon bilinear operators A\ and Ni are rewritten as 

A\ = Ui{q)vi{q) + u 2 i (q)Al{q) - v 2 i {q)A i {q) - Ui{q) Vi {q)N i(q), 
Ni = 2v 2 (q) + 2u l {q)v l {q){A\{q) + A t (q)) + K^) " v 2 (q))Ni(q). 



(4-9) 
(4-10) 



(4-11) 
(4-12) 



The quasiparticle bilinear operators A\ (q) , Ai (q) , and iVj (q) satisfy the following 
commutation relations: 



Ai{q),A\,{q)\ =8 iV {\-N i (q)), 
Ni(q),A\,(qj\ =25 u ,A\,(q). 



(4-13) 
(4-14) 



The particle number No, the quadrupole deformation D(q), the monopole-pairing 
gap Z^o (q) , and the quadrupole-pairing gap A 2 (q) are given by the expectation values 
with respect to the mean-field state vector \4>(q)): 



N = (ct>(q)\N\ct>(q))=2j2vI(q), 

i>0 

D(q) = (<f>(q)\ D = 2Y,d t a l v 2 (q), 

i>0 

A (q) = Go {<f>(q)\ i f \<j>(q)) = Go^^iQMQ), 

i>0 

A 2 (q) = G 2 (m\B ] \m) = G 2 Y, diaMqMq). 



(4-15) 
(4-16) 
(4-17) 
(4-18) 



i>0 



Below, we often omit the (/-dependence in expressions, for example, writing 
Ai(q) as Ai. It should be kept in mind, however, that these quantities are locally 
defined with respect to the quasiparticle vacuum \(/>(q)) and depend on q. 



4.3. The ASCC equations for separable interactions 

The ASCC equations for the separable Hamiltonian are given by 

5{<f>(q)\h M (q)\ct>(q))=0, 



(4-19) 



HHq)\ [h M (q),Q(q)} - E4^i _) - l - B ^)P(q) \m) = 0, (4-20) 



HHq)\ 



h M (q),-B(q)P(q) 



E fpJ(Q)^ +) - B(q)C(q)Q(q) ~ £ /g!(<z)£ (+) 



(+)/ 
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\[[h M {q)^Q{q)]M)] 



f, 



(-) 



f N (q)N\0(q))=O, 



(4-21) 



where /im(<?) denotes the self-consistent mean-field Hamiltonian in the moving frame, 
defined by 

h M (q) = h(q) - X(q)N - ^Q(?), (4-22) 

with 



h(q) = h -^ Ks FM F« |0( g )). 



We also define the following quantities 

&9) = -n. (<Kq)\[F$-\Q(q)]\<Kq)), 
f£(q) = k s (0(g)| [F S M ±B(q)P(q)] |0(g)) 



IP,s 



2^ 



/jv(g) = B(q) 



dX 
dq ' 



Q(q) 



Note that all matrix elements are real, so that (4>(q) \ Fs \4>{q)) = 0. The fifth term 
in Eq. (|4-2ip appears only in the gauge-invariant formulation of the ASCC equations 
because there is contribution from this term only if the Q(q) contains the -B-part. 

4.4. The moving-frame HFB equation 

The moving-frame HFB equation (I4-19D at a given q determines the time-even 
TDHB state vector \cp(q)). The variation in Eq. (|4-19p is taken with respect to 
arbitrary two quasiparticle states: 

5 |0(g)) =aj(g)at(g) |0(g)). (4-28) 

If we know the operator Q(q), we can solve this equation using the gradient method 
so as to eliminate the two-quasiparticle terms proportional to A\ and Aj. The 
quantities, A(g) and dV/dq, can be regarded as Lagrange multipliers which are de- 
termined by the following two constraints. The first is the particle number constraint 
given by (14-15p . This constraint specifies the location in the particle number space. 
The second constraint is written as f)3- 14|) . For the Q operator defined by (|3-10p . 
this equation yields 

(<P(q)\Q(q - Sq) \<f>(q)) = 2^Qi(<Z - 6q)(vi(q) 2 - Vi (q - 5q) 2 ) = 5q (4-29) 

j>0 

4.5. The local harmonic equations 

We solve the local harmonic equations to obtain operators Q(q) and P(q)- The 
collective coordinate operator Q{q) is written as 

Q(q) = £ Gifa) : & ■= E {Qt(q)(A] + A t ) + Qf(q)Ni] , (4-30) 

i i>0 



\Hq)) , 



(4-23) 

(4-24) 
(4-25) 

(4-26) 
(4-27) 
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while the collective momentum operator P(q) is expressed as 

P{q)=iY,P,{q)(A\-A l ). (4-31) 

As mentioned in the preceding section, the S-part of the operator P(q) is unneces- 
sary in the second order with respect to the collective momentum p. 

We solve the local harmonic equations in the following way: Assume that the 
solution Q^ n ~ l \q), obtained in the previous iteration step, of the local harmonic 
equations and the solutions, such as \4>{q)) and V(q), of the moving-frame HFB 
equation are available. (The superscript n is omitted except for Q(q).) In solving 
the local harmonic equations, we note that the moving-frame Hamiltonian /im(<?) 
and the operators F^ are expressed in terms of the quasiparticle bilinear operators 
A\ , Ai , and iV j as 

h M (q) = V(q) + Y,E l (q)N u (4-32) 

Fj+) = Mq)\ jfK+> |<K<Z)> + f$ + fj£ 

= \<f>(q)) + ^Fi+ ) (,)(Aj + Ai) + ^4+ ) (,)iV l , (4-33) 

i>0 i>0 

Fs ( - ) =J2 F As(. i )( A l-A i ). (4-34) 



i>0 



Here, 



^2(0 = ^-*?), F$(i) = ±d i a i (uf-v 2 ), F^(i) = 2d i a l u l v l , (4-35) 



^l~l ) (0 = ~, F i J(i) = ~d l a l , F { J(i) = 0, (4-36) 

F B^l(^ = ~ u i v i-, ^bI^S) = -diCFiUiVi, F^l (i) = di(Ji(u 2 - v 2 ), 

(4-37) 



Ei{q) = (u 2 - v 2 ) (a - xdmD{q) - \{q) - l) \q) \ - 2{A (q) + d^ 2 (g)) W 

(4-38) 

These quantities are determined by solving the moving- frame HFB equation (|4- 19[) . 
For later convenience, we define the following quasiparticle bilinear operators: 

= [F^\^Q^- 1 Hq)] = 2j2R^ s (i)(Al T A i ), (4-39) 

with 

<]» = f «(0<tf '"""to - F^i)Qf^%)) . (4-40) 

<]» = -^^ ) BQf < - 1, W (4-4D 
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We can express the matrix elements and P{ in terms of /qJ, /pV? /r~ s an d 

/at by substituting Eqs. ([¥30]) and (|¥3T]) into Eqs. (gEHD and (14^211) : 



1 



+ 



dq 

2E, 



(2^) 2 -^ 2 + 2^gf (n ^ 1) 

(4-42) 



9g 

^^-J^qf.^ E^'Wifc'. (-3) 

where 

iVi = 2u i ( 9 )i; i (g), (4-44) 

/^. = /ft ) (g) + /fc ) (g). ( 4 - 45 ) 
o> = y/B{q)C{q). (4-46) 

Substituting Eqs. (JS3QD, (l¥3TD and (g^g) into Eqs. (ITOj) . (I¥25l) and g26|), we 
obtain 

4i=2^E^W (n) , ( 4 - 47 ) 

= E {i^Wfl + W } • (4-48) 

Note that /q 3 = 0. From the canonical variable condition, the orthogonality of the 
collective and number fluctuation modes is required: 

^(q)\[N,P(q)Mq)) = 2iJ2NiPi = 0. (4-49) 

i>0 

Eliminating Qf (n) and Pi from Eqs. (033), (H3HD, and fljggp with use of Eqs. (g32j) 
and ()4-43p . we finally obtain the dispersion equation 

S(^ 2 )-/=0, (4-50) 

for the quantity /= f{q) = {fq~}, fg^, fpR >v fpR >2 , f { pR i3 , In}- Here S = {,%} is 
a 6 x 6 matrix whose elements are given by 



u? - 2%E l Q B( - n ~ 1) 





= 2« 1 {^)(J^,FH) 


-2S^(F^,R^)}-1, 


(4-51a) 


S12 


= 2 Kl {s^(F^,F^) 


-2S^(F^,R^)}, 


(4-51b) 


Sl 3 


= 2K lS W(Fti,F$), 




(4-51c) 


S\4 


= 2 K1 SW(F£IF$), 




(4-51d) 


S\5 


= 2 K1 SW(F£IF$), 




(4-51e) 


S\6 


= 2 Kl SW(F { -lN), 




(4-51f) 
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o 

<J21 




-2S^ >{F\J,R\{)j, 


(A rn„\ 

(4-oza) 


£22 


= 2 K2 {S^(F^IF^) 




(4-52b) 


-523 


= 2 K2 SW(F$,F$), 




(4-52c) 


#24 


= 2n 2 S^(F^,F^), 




(4-52d) 


S25 


= 2n 2 S^\F^,F^), 




(4-52e) 


S26 


= 2 K2 S^(F^,N), 




(4-52f) 



531 = 2/ci 

532 = 2/Ci 

#33 = 2/Ci 
S34 = 2/Ci 
£35 = 2/Ci 

Ssq = 2ki 

541 = 2k 2 

542 = 2k 2 

543 = 2k 2 

544 = 2k 2 

545 = 2^2 
S4Q = 2k 2 

#51 = 2^3 

552 = 2k 3 

553 = 2«3 



! S^(F^,F^)-S 



S {1) (F$,F$) + S W 
SM(F$,F$) + S V 
SV(F$,F$) + S W 
S^(F^,N) + S^(R 



u^(F^,F^)-S 



u 2 S^ 2 \F^,F^ 



s 



s (1 \f£If a1 



A,2 ' A A,2 

<{+)\ +S (2) 



SM(F$,F$) + SW 
SV(F$,F$) + SW 
SW{F$,N) + SW(R 



^\f$,f£)-s 

u^\f£If£)-S 
SM(F%,F$) + S W 



HF A :t,Rti) + sU(R { 2lF£ 



-h 

2 , 



1} (^2 ' fl i,2 ) + 5 ' (1) ' F i,2 ) 



o(+) 

a A,2' r A,l 

o(+) p (+) 
^^4,2' ^A.2 



i? 



(+) P (+) 



7?^ 

A,2' A,3 



(+) p (+)n\ 

A,3' r A,l ;/ ' 



i? 



-2^)(<;,4;)) 

(4-53a) 

- 2 ^)(<;,4;)) 

(4-53b) 
(4-53c) 

(4-53d) 

(4-53e) 

(4-53f) 

-2S< 2 >(*« 
(4-54a) 

- 2S ( - 2 \r ( -^1, R Ai l) 
(4- 54b) 

(4-54c) 
(4-54d) 
(4-54e) 
(4-54f) 

-2^)(<),<)) 
(4-55a) 

-2^)(<),<)) 
(4-55b) 
(4-55c) 
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5 54 = 2k 3 {s (1 \F^If^)+S^\R%F^)} , (4-55d) 

5 55 = 2k 3 {S^HF^IF^)+S( 2 \R%F^)} - 1, (4-55e) 

5 56 = 2 K3 {SW(F$,N) + 5( 2 )(i?«,iV)} , (4-55f) 



5 6 i = u?S<- 2 \N,F$) - S^(N,R { -}), S 62 = u 2 S&(N,F$) - sW(N,R$), 

(4-56a) 

5 63 = S« (AT, F« ) , S 64 = 5« (iV, ) , (4- 56b) 

5 65 = 5«(iV,FW), S m = S^(N,N). (4-56c) 

Here, the quantities and S^ 2 ) are defined by 

SM(X,Y) = y ——XiYi, (4-57) 

(X, Y ) = V ,„ XiYi. (4-58) 

t (2E i ( Q )y--H<l) + 2^E l Q\ B{n - 1)) 

The unknown quantities in the dispersion equation (|4-50p are f(q) and uj 2 (q). The 
squared frequency uj 2 (q) can be determined by the condition that the matrix S(uj 2 (q)) 
has no inverse: 

det S{u 2 (q)) = 0. (4-59) 

In the case that there are many solutions u 2 (q) satisfying this equation, we choose 
the smallest of these (including negative values) as the collective mode. Once the 
value of oj 2 {q) and, consequently, the matrix S(q) is specified, the direction of the 
vector f(q) is found. Then, its absolute value is fixed by the normalization condition 
for the collective mode, i.e., 



(HQ)\ [Q {n) (q),P(q)] \<t>(q)) = Q?™ WW = L ^ 

i>0 



The choice of the signs of (q) and P(q) are still arbitrary. This sign specifies 
the "rear" and "front" of the one-dimensional collective path. 

The i?-part of Q{q) is automatically determined in terms of its A-part according 
to its definition (g^J): 

Qf {n \q) = ^^Qt {n \q). (4-61) 

4.6. Gauge fixing 

Under the gauge transformation f)3- 5[) , the quantities, fk I (q) , fk \ (q) , and fjy (q) 
appearing in the dispersion equation (|4-50p transform as 

(4-62) 

f { ~ ] 



^qM^-^Mq), (4-63) 
f N (q) ~ au; 2 (q). (4-64) 
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These properties clearly indicate that one of the above three quantities can be elim- 
inated: By choosing an appropriate value for a (gauge fixing), we can reduce the 
dimension of the dispersion equation (|4-50p to a 5 x 5 matrix equation. Namely, 
Eq. (|4-50p is redundant and the gauge fixing is equivalent to the reduction of its 
dimension. The QRPA gauge corresponds to setting /n(q) = 0. In our previous pa- 
per 57 ) we set fg I (q) = which corresponds to another gauge. Because the quantity 

/q i(q) represents the contribution from the time-odd component of the monopole- 
pairing interaction, let us call this gauge "ETOP (eliminating time-odd pairing) 
gauge" 

4.7. Requantization 

The solution of the ASCC equations yields the classical collective Hamiltonian: 



We then obtain the quantum collective Hamiltonian by carrying out the canonical 



quantization Tt(q,p) — > 7i I q, — -ft I • Note that, in this quantization step, there is 



no ambiguity associated with the ordering of q and p, because the coordinate scale 
is chosen such that the inverse mass function is unity, i.e., B(q) = 1. 

§5. Numerical test of internal consistency of the proposed scheme 

5.1. Details of numerical calculation 

We solve the gauge- invariant ASCC equations for the multi-0(4) model with the 
same parameters as in the previous papers. 57 - 1 ' 62 ) The system consists of 28 particles 
(one kind of Fermion). The model space consists of three j-shells, labeled ji, j'2, J3, 
with pair degeneracies f2j i = 14, Qj 2 = 10, f2j 3 = 4, single-particle energies = 
0, ej 2 = 1.0, ej 3 = 3.5, and the single-particle quadrupole moments dj 1 = 2, dj 2 = 
dj 3 = 1. Within this model space, the deformation D = (<f)(q) \ D \4>(q)) ranges from 
D min = -42 to D max — 42. The calculation is done for the quadrupole-interaction 
strength x = 0-04 and the monopole-pairing-interaction strengths Go = 0.14,0.16 
and 0.20. The properties of the system change from the double-well situation (Go = 
0.14) to the spherical vibrator (Go = 0.20) according to the value of Go- The effect 
of the quadrupole pairing is studied by comparing the results for G2 = 0.00, 0.02 
and 0.04. As pointed out in Ref. 57) the quadrupole pairing gives significant effects 
on the collective mass. It is not essential, however, for the discussion on the gauge- 
fixing condition. Therefore, we present the results for G2 = in the next subsection 
and show its effect in the final subsection only. The calculation starts from one 
of the HFB equilibrium state labeled by q = (see Fig. [2D. For the deformed 
cases (Go = 0.14 and 0.16), the HFB equilibrium state having positive (prolate) 
deformation is chosen as a starting point. 

5.2. Comparison of the two gauge fixing conditions 

Let us examine whether or not we can find the gauge independent solution of 
the ASCC equations. Existence of the collective path that simultaneously satisfies 
all equations of the ASCC method is not self-evident. The aim of the numerical cal- 
culation here is to check internal consistency of the equations set up in the preceding 
section. We solve the gauge-invariant ASCC equations with two different gauge fix- 



n(q,p) = -p 2 + V(q). 



(4-65) 
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ing conditions: the QRPA gauge (/jv(<7) = 0) and the ETOP gauge (fqliq) = 0). 
In the QRPA gauge, the chemical potential X(q) along the collective path is set to 
be constant, while in the ETOP gauge, the time-odd contribution of the monopole 
pairing interaction is fully eliminated from the ASCC equations. 

Figure [1] shows the collective potential V(q) and the monopole pairing gap Ao(q) 
as functions of the quadrupole deformation D(q). Figure [5] displays the relation 
between the collective variables q and the quadrupole deformation D as well as the 
squared frequency oj 2 {q) obtained by solving the local-harmonic equations. The 
collective mass M(D(q)) = (dq/dDf is plotted as a function of D in Fig. [3j We 
find that the calculation using the ETOP gauge works very well and the collective 
path connecting the two local (oblate and prolate) minima with different signs of 
deformation are successfully obtained. In contrast, the calculation using the QRPA 
gauge encounters a point where we cannot proceed any more. In the region where 
the solutions have been found for both gauges, they well agree with each other. It 
should be the case because these are gauge invariant quantities. The cause of the 
difficulty encountered in the QRPA gauge is understood in the following analysis. 

Figures H] and [5] display the chemical potential X(q) and the quantities, fk l(q) 
and /at((?), respectively. Their values depend on the gauge adopted. If the QRPA 
gauge is used, X(q) should be constant along the path because of the condition 
fN(q) = B(q)d\/dq = 0. We find, however, X(q) diverges near the inflection point of 
the collective potential where u(q) 2 = d 2 V/dq 2 = 0. This divergence occurs because 
the inflection point is a singularity for the gauge transformation (|4-64p where an 
arbitrary a gives the same /at(q). Thus, the calculation using the QRPA gauge stops 
at the inflection point. On the other hand, we can go over the inflection point using 
the ETOP gauge, because the gauge transformation for fk ^(g) (|4-62[) involves only 
the monopole pairing gap Ao(q) which always takes finite values along the collective 
path (except at the limit of the model space). In these figures we also present the 
results that are obtained by the following procedure: After determining the collective 
paths with the use of the ETOP gauge, we calculate the gauge-dependent quantities, 
•M<?)> Iq i (<?)> /jv(?)) by switching to the QRPA gauge using the relations 

.(-) _ 4/jy(g)(ETOP)A)(g) 

/Q,i(QRPA) - w 2( g ) l b ' la J 



^(<?)(QRPA) =A(<?)(ETOP) -7l7Z\ — ( 5 ' lb ) 



4/at(<?)(etop) dV 
0J 2 {q) dq 

We see in Figs. 0] and [5] that the results obtained by this procedure agree with those 
calculated by using the QRPA gauge (/at(<?) = 0) from the beginning (in the region 
of deformation D(q) where the collective path can be obtained using the QRPA 
gauge). This agreement demonstrates that the collective paths determined by using 
different gauge fixing conditions are the same, as it should be. Nevertheless, there 
is a suitable gauge fixing condition in finding solutions of the ASCC equations and 
constructing the collective path. For the multi-0(4) model with superfluidity, we 
find that the ETOP gauge is more useful than the QRPA gauge, because the gauge 
transformation (|4-62p is well defined as long as the pairing gap Ao(q) is non-zero. 

5.3. Comparison with the previous calculation 

In our previous paper, 57 ) we employed the ETOP gauge condition (fqliq) = 
0) but the i?-part of Q(q) is ignored. Let us evaluate the error caused by this 
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approximation. The results of such calculation are presented also in Figs. [T][4] and 
compared with those of the full calculation. We see that they agree well indicating 
that the approximation of ignoring the i3-part is rather good. 
In Fig. [5l we present the quantity 

f { Q,!(l) = (<Ko)\ i^'\Q(Q)} l<K<?)> = (5-2) 

i>0 

evaluated using the Q(q) operator that is obtained by ignoring the -B-part in the 
process of solving the ASCC equations. This quantity should be zero if the Q(q) 
operator determined by the gauge-invariant ASCC equations is used. We see that 
the deviation from zero is negligible (except near the limit of the model space) , again 
indicating that the approximation is good. 

The quantum spectra and transition strengths are displayed in Fig. [6J These 
are obtained by solving the Schrodinger equation for the quantized collective Hamil- 
tonian. In this figure, effects of the quadrupole pairing interaction are also shown. 
We see that the results of the previous calculation (in which the -B-part of Q(q) is 
ignored) are quite similar to those of full calculation (including the B-part), and 
both results well reproduce the trend of the excitation spectra obtained by exact 
diagonalization of the microscopic Hamiltonian (Fig. E]). The numerical calculation 
presented above thus suggest that the approximation of ignoring the £?-part of Q(q), 
adopted in our previous paper, 57 ) is justified and it may serve as an economical way 
of determining the collective path. 

§6. Concluding Remarks 

We have shown that the basic equations of the ASCC method are invariant 
against transformations involving the angle in the gauge space conjugate to the 
particle-number. By virtue of this invariance, a clean separation between the large- 
amplitude collective motion and the pairing rotational motion can be achieved, en- 
abling us to restore the particle-number symmetry broken by the HFB approxima- 
tion. We have formulated the ASCC method explicitly in a gauge-invariant form, 
then, applied it to the multi-0(4) model using different gauge-fixing procedures. The 
calculations using different gauges indeed yield the same results for gauge-invariant 
quantities, such as the collective path, the collective mass parameter, and spectra 
obtained by the requantization of the collective Hamiltonian. We have suggested a 
gauge-fixing prescription that can be used in realistic calculations. 

The explicit gauge invariance requires a 5-part (a^a part) of the collective coor- 
dinate operator Q(q). Actually, Ref. 27) remarks that the separation of the Nambu- 
Goldstone modes in the local harmonic formulation requires higher-order terms in 
the collective coordinate. This is consistent with the present conclusion in the gauge- 
invariant formalism. We have also demonstrated that the approximation to neglect 
the B-part leads to results almost identical to those of the full calculation, at least 
for the multi-0(4) model. 

We are now investigating the oblate-prolate shape coexistence phenomena 54 ) -56 ) 
in nuclei around 68 Se using the pairing-plus-quadrupole interactions 63 ) -65 ) using the 
prescription on the basis of the new formulation of the ASCC method proposed in 
this paper. The result will be reported in a near future. 
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. 1. Collective potentials V(q) and monopole pairing gaps Ao(q) plotted as functions of the 
quadrupole deformation D. The upper, middle and lower panels display the results for Go = 
0.14,0.16 and 0.20, respectively. In each column, results of different calculation are compared; 
those obtaind using the QRPA gauge (/jv(<?) = 0) and the ETOP gauge (/q7i( ( ?) ~ 0) are 
plotted by solid (red) and dotted (blue) lines, respectively, while those obtaind ignoring the 
B-part of Q(q) (i.e., putting Qf (q) = 0) by dashed lines. 
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Fig. 2. Left column: Relation between the collective coordinate q and the quadrupole deformation 
D(q) = {4>(q)\ D \(f>(q)) . The point q — corresponds to the HFB equilibrium, which is the 
starting point of the numerical calculation. Right column: Squared frequencies u) 2 (q) of the 
local harmonic equations, plotted as functions of D. Note that they are negative; i.e., u>(q) is 
purely imaginary in the region where the curvature of the collective potential is negative. The 
upper, middle and lower rows display the results for Go = 0.14, 0.16 and 0.20, respectively. See 
the caption of Fig. [T] 
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. 3. The ASCC collective mass M(q(D)) as functions of the deformation D. The upper, middle 
and lower rows display the results for Go = 0.14, 0.16 and 0.20, respectively. See the caption 
of Fig. ED 
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Fig. 4. The chemical potentials X(q) as functions of the deformation D. The upper, middle and 
lower rows display the results for Go = 0.14, 0.16 and 0.20, respectively. In each column, results 
of different calculation are compared; The solid (red) and dotted (blue) lines indicate the results 
obtaind using the QRPA gauge (/iv(g) = 0) and the ETOP gauge (fo~i(q) = 0), respectively, 
while the dashed lines indicate those obtaind using the ETOP gauge but ignoring the B-part of 
Q(q) (i.e., putting Qf(q) = 0). The dotted (purple) lines labeled "combination" indicate the 
results calculated by switching to the QRPA gauge after the collective path is determined by 
using the ETOP gauge. 
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functions of the deformation D. The upper, middle and lower rows display the results for 
Go = 0.14, 0.16 and 0.20, respectively. See the caption of Fig. [4] 
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Fig. 6. Excitation spectra and quadrupole transition matrix elements, left column: The result 
obtaind using the ETOP gauge (fg I (q) = 0) but ignoring the B-part of the collective coordinate 
operator Q(q). This is the same as presented in Ref. 57). right column: The result obtained by 
solving the gauge-invariant ASCC equations using the ETOP gauge. The upper, middle and 
lower rows display the results for Go = 0.14, 0.16 and 0.20, respectively. In each row, the results 
for G2 = 0.00, 0.02 and 0.04 are compared. The numbers adjacent to the vertical lines are the 
absolute values of the transition matrix elements. The matrix elements between the doublets 
are indicated beside them. 
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Fig. 7. Excitation spectra calculated with the exact diagonalization of the multi-0(4) model Hamil- 
tonian. The upper, middle and lower rows display the results for Go = 0.14, 0.16, and 0.20, 
respectively. In each row, the results for G2 = 0.00, 0.02, 0.04 are compared. The numbers 
adjacent to vertical lines indicate the transition matrix elements. The matrix elements between 
the doublets are indicated beside them. 



